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ABSTRACT 

A key parameter in the design trade-offs made during AlGaN/GaN HEMTs development for microwave power 
amplifiers is the channel temperature. An accurate determination can, in general, only be found using detailed 
software; however, a quick estimate is always helpful, as it speeds up the design cycle. This paper gives a simple 
technique to estimate the channel temperature of a generic microwave AlGaN/GaN HEMT on SiC or Sapphire, 
while incorporating the temperature dependence of the thermal conductivity. The procedure is validated by 
comparing its predictions with the experimentally measured temperatures in microwave devices presented in three 
recently published articles. The model predicts the temperature to within 5 to 10 percent of the true average channel 
temperature. The calculation strategy is extended to determine device temperature in power combining MMICs for 
solid-state power amplifiers (SSPAs). 


I. INTRODUCTION 

At the present time many research groups are studying AlGaN/GaN HEMTs for high frequency (lip to 
Ka-band), high power microwave applications [1-4]. The large RF power density generated in these amplifiers 
causes considerable self-heating, and a reasonably accurate estimate of the channel temperature is often desired. 
Knowledge of the temperature is essential as both the carrier mobility and reliability suffer. The mobility reduces 
with increasing temperature as (1/T) 2 \ with resulting decrease in DC and RF performance [5]. Both short and long 
term reliability are key parameters, which also reflect the need to rapidly calculate the temperature profile in a 
particular design. 

Both the long and short period output power changes measured in many devices under stress testing are related 
to the device temperature. The trade-off in the number of fingers in a MMIC cell is controlled to some extent by the 
thermal hot spots that develop near the edge of the gate in the gate-drain access regions in the center of a cell. 
Several 2- and 3-D simulators [6-9] are available to determine temperature profiles. A 2-D electro-thermal package, 
MINIMOS-NT [10] considers the coupling between the current and heat transport equations, and thereby performs a 
self-consistent calculation. The SILVACO [11] suite of software also performs self-consistent calculations. A few 
recent articles [12-14] give techniques to quickly ascertain junction temperatures under certain conditions; two of 
which neglect the temperature dependence of the thermal conductivity. This article demonstrates a procedure to 
estimate the channel temperature in microwave AlGaN/GaN HEMTs on either SiC or Sapphire, using simple 
analytic expressions, which includes the change in thermal conductivity with temperature. The device is modeled 
using the thermal resistance method, and is applicable to multi-gate, multi-material, structures, and a wide range of 
power dissipation levels. 

We do not want to mislead the reader into thinking we have a nearly final answer to the many problems that 
plague AlGaN/GaN HEMTs (such as current collapse, permanent electrical change after initial electrothermal stress, 
etc.). What we demonstrate is the reasonably accurate temperature estimate by using the quiescent dissipated power. 
This is a contribution in lieu of the extreme difficulty in modeling the device accurately using existing software. As 
stated in [13] the electrothermal modeling of devices is such a difficult problem that the present models (circa 2001) 
are either oversimplified or rather inefficient from a computational point of view. For example, Wachutka [15] 
developed the first rigorous specification of heat generation in a semiconductor based on thermodynamic concepts. 
This model is available with the SILVACO [11] package, whereas the simpler expression, H = J E is also 
available to reduce complexity and simulation time. Efficient thermal management of electronic components is 
essential to minimize the influence of thermomechanically induced stress and thermal load. The fact that the 
2D-electron gas in AlGaN/GaN HEMTs is abut 50 percent supported by piezoelectric polarization due to 
mechanical strain; the interaction of temperature, strain, piezoelectric fields, and carrier densities at critical positions 
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in the device is of paramount importance. Here we take a less exact approach to demonstrate the application of 
simple formulas for the analysis of combined devices used for solid state power amplifiers (SSPAs). We do not 
attempt to ascertain the complex interaction in a power device, just the average temperature to determine reasonable 
spacing between electrically combined devices as well as the package layout. 


II. ANALYSIS PROCEDURE 

A generic microwave power AlGaN/GaN HEMT layout and its cross section are shown in figures 1 and 2, 
respectively. Being a power device it has a multi-finger gate structure, and the die is attached to a Cu-tungsten heat 
sink for efficient thermal management. The thermal conductivities for the many material layers in the cross-section 
were determined by making rough averages of the values reported in the literature. To convert the actual device 
depicted in figures 1 and 2, into one wherein analytical formulas are available, we use the following approach. On 
the left in figure 3 we have shown an equivalent effective heat aperture to model the heat generated on the top 
surface (primarily in the gate-drain access regions). Then on the right we show a series of thermal resistances which 
depend on the thermal conductivity and thickness of each layer of the device. In cases where the separation between 
the gate fingers is large with respect to the chip thickness, the heat source at the position we show as /channel in the 
figure, is a thin strip of uniform heat generation, which is about the same physical size as the gate electrode, and 
located in the gate -drain access region. The technique consists of determining the thermal footprints of the heat 
sources at each successive lower layer, and then determining that layer’s thermal resistance. 



Cu-tungsten carrier plate 


Figure 1 . — Schematic of an 8 gate AlGaN/GaN HEMT 
for a microwave power amplifer. 
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Figure 2. — Cross-section model of the HEMT device, 
showing approximate layers of material with aver- 
age, or ranges, of thermal conductivities. 
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Effective heat 
aperture -a 




Figure 3. — (a) Equivalent heat source of the 8 
gate HEMT device for modeling thermal profile, 
(b) Equivalent circuit model representation of 
the thermal resistances of the HEMT device. 


A. Expressions for Temperature Dependent Thermal Conductivity 
In principle, in any heat transport problem, one solves for the temperature field from the heat flow equation 

V '{ K ; (T)vr} = -P (1) 

where k ,( 7) is the temperature dependent thermal conductivity in the i-th layer, T is the temperature in Kelvin, and 
P is the dissipated power density in W/m 3 . The equation is altered by defining an equivalent artificial temperature 
via the Kirchhoff transformation 

(2) 

where T„ is the reference temperature in K, Tf is the transformed temperature in the i-th layer, and T is the real 
temperature there in K. The transform linearizes equation (1) with the result 

K/(7’ 0 )V 2 7} e =-P (3) 

While the Kirchhoff transform linearizes the heat equation, the nonlinearities inherent in equation (1) become a 
problem when layers of different materials are considered. When the layers have different functional relationships 
for their thermal conductivities k,(7), then the transformed temperature across a boundary separating two 
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homogeneous regions will be discontinuous. Only when the ratio of the conductivities k,(T)/k !+ i(7) is temperature 
independent, will the transformed temperature be continuous across the interface [16]. With this in mind, our first 
task is to determine analytical expressions for the temperature dependent thermal conductivities of the GaN, SiC, 
and Sapphire layers. We have not found an analytical expression for that of GaN which was developed from 
measured data, so we choose the following expression 


K GaN ( 7 ’)- 1 -6 



W 

cm-K 


( 4 ) 


This was arrived at by assuming GaN is sufficiently similar to GaP [17-18] (to choose the exponent), and the 
coefficient is in the range reported in the literature. For SiC we choose [19] 



W 

cm-K 


( 5 ) 


which was obtained from curve fits to experimental data. For Sapphire we assume [20] 

K( r),^_ J 04<, >* 300 ) w <6, 

v ’ T - 159 T cm-K 

The first part of expression (6) was determined from curve fits of measured data, and the second one was an attempt 
to have a 1/T variation. Figure 4 plots the above thermal conductivities. Observation of equations (4) to (6) shows 
that the ratios of the thermal conductivities are not independent of temperature, but they are not real strong functions 
of temperature. The change in the ratio for GaN on SiC is 5 percent over the temperature range of 300 to 500K. The 
corresponding change for GaN on Sapphire is 23 percent. This variability is not a problem as we determine the 
actual temperature at each interface at each step in the calculation. The variability becomes a problem when 
developing computational codes. We recognize that the leading coefficients and exponents are approximate as they 
are dependent on varying material parameters such as percent of cubic (Zincblende) phase in an otherwise hex 
(Wurtzite) phase sample. The dislocation densities, cracks, and other growth dependent factors will force changes in 
the formulas, but they are accurate enough for practical purposes. 



Figure 4. — The change in thermal conductivity 
versus temperature for SiC and sapphire 
substrates, and GaN epitaxial layers. 
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B. Expressions for Thermal Resistance for Single Heat Source (Gate Finger) 

The thermal resistance R TH in (K/W) is determined by using the complete 3-D nature of the heat spreading. 
A device with a single gate finger can be modeled as a single heat source of length 21 x and width 2/,, as shown in 
figure 5. The thermal resistance of any layer in the stack is given by [12]. 


where 


R m ~ ~ 
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Figure 5. — Geometry for calculating thermal 
resistance for (a) a device with a single 
gate finger, (b) multi-gate structures. 
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where 



(9a) 

(9b) 

(9c) 

(9d) 

(9e) 


where the angles a and (3, the lengths L x and L v , and the thickness w are as shown in figure 5(a). This expression for 
Rtii is a general one, applicable to any layer in the resistance stack in figure 2. 


C. Expression for Thermal Resistance for Multiple Heat Sources (Multi-Gate Structures) 

In a multi-gate structure there is a strong possibility that the fluxes will intersect. The thermal resistance of any 
layer in the stack is given by [21] (with typos corrected) 


where 

n the number of fingers 
Z total gate width = n(2l y ) 


R m - 


n 

Zkti 

2(«-l) («-2) 

InM In P 


and 



( 10 ) 


( 11 ) 


( 12 ) 


The step-by-step procedure is as follows. Start with the heat sources on the top surface; which are stripes with 
the short side dimension being that of the gate length, or the gate-drain spacing. Using equations (8) and (9), 
determine the flow pattern as shown in figure 5(a). Use the room temperature values of k, as the correct values to 
use are not yet available, since the temperature on the top surface is unknown at this point. Then the footprint on the 
next lower surface is determined. Continuing through all layers, one determines the thermal resistance stack. The 
calculation continues by starting at the assumed Cu-tungsten heat sink temperature and working upward through the 
layers. The change in temperature across a layer is given by 
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AT i = R TH i P DISS 


(13) 


where P D1SS is the dissipated power in watts, i.e. the I D V D product, AT] = the temperature rise in K, and R TH = the 

thermal resistance of the layer. The actual temperatures are used in the layers with constant thermal conductivity, 
and the artificial ones are determined in the GaN, SiC, and Sapphire layers. The actual temperature T A at the surface 
of these materials is found from inverting the Kirchhoff transform, which yields 


p A (GaN) = 



^(SiC) = 



T A (Sapphire) = T 0 exp 



(14) 

(15) 

(16) 


where T a is the reference temperature, (the actual temperature at the lower surface of a given layer). With the first 
estimate of the channel temperature available, one may iterate using a more accurate value for k(T), rather than the 
room temperature value assumed. However, we find that this error reduction step is not usually needed. The 
procedure is now illustrated in detail in section 111(A). 


III. COMPARISON WITH EXPERIMENTAL DATA 
A. Single Gate Device on SiC 



Figure 6. — Schematic of an equivalent heat source 
for a single gate device used in the measurement 
of the channel temperature by Raman spectro- 
scopy [22]. 


For a device with a single gate of length 2 l x = 1 pm 
and width 2 l y = 200 pm, we model the heat source as a 
1x200 pm stripe as shown in figure 6. We neglect the 
0.28 pm AlGaN layer, because it is thin, and in [22] this 
was neglected in the calculation. 

First steps: Calculate the thermal footprints at all 
material interfaces using equations (8) and (9). The 
1x200 pm source region on the top of the GaN becomes 
a 3.1x200.2 pm source on the GaN/SiC interface. Using 
this we find that the footprint at the SiC/heatsink plane 
spreads to 673 pm (which is greater than the 72 pm 
dimension of the chip). Thus edge effects come into play, 
and may be handled using the techniques in [24]. We 
however, neglect these corrections, which we assume 
will be small in this particular case. 

Next we determine the thermal resistance of the 
layers. Using equations (7) to (9) we find the first value for 
R m (GaN) to be 3. 14 K/W. We also use equation (10) with 
n = 1, and obtain RthZk = 0.587. Assume k(T) =1.6 (again 
the 300K value which we used in equations (7) to (9)) as 
we don’t know the actual temperature and thus k for that 
temperature. This will very slightly underestimate the final 
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channel temperature. We find R m (GaN) = 18.35 K/W. Notice that this value is larger than the first value by a factor 
of 5.86. We assume this value is less reliable, and use the previously calculated one, i.e. R m = 3.14 K/W. Using 
equations (7) and (8) we find the footprint on the SiC layer is 3.1x200.2 pm. Use equation (10) again for this layer 
and obtain R m (SiC) = 29.72 K/W. The dissipated power is 1.754 W which yields AT(SiC) = 52 K, AT(GaN) = 

32.2 K. Using equation (15) gives T(SiC, actual) = 86.6 °C. Then in the GaN layer, T = 359.6 + 32.2 = 391.8 K, 
which is the modeled channel temperature, Tchannel (GaN, actual) = 121 °C. In [22] and [23], temperature 
measurements on microwave AlGaN/GaN HEMTs on both SiC or sapphire were made using Raman spectroscopy. 
The spatial resolution was 1 pm and the temperature accuracy was within 10 °C. The measured value for the channel 
temperature Tchannel was 124±5 °C. The 3-D model used in [22] assumed a 1/T variation in k and predicted 
Tchannel = 99 °C. The reason for the low accuracy is perhaps due to the assumed variation of k with temperature. 
Thus our model achieves much better results than the 3-D model. 


B. Single Gate Device on Sapphire 

For the sapphire substrate case we choose 2 l x = 4 pm (gate length), k (sapphire, 300 K) = 0.49; then calculate 
R m (GaN) = 8.35 K/W, R m (SiC) = 180.5 K/W, AT(sapphire) = 1 17.3 K, and AT(GaN) = 5.4 K. Then the calculated 
temperature T CH annel using our model on the top surface becomes 176 °C. The measured value was 180 °C and the 
3-D model predicted 140 °C. Hence our method again predicts the temperature more accurately. 


C. Multi-Finger GaN Device on SiC 

An 8-gate device on SiC (fig. 1) for 2 to 4 GHz amplifier applications was measured in [23], and the thermal 
profile is shown in figure 7. Using the outlined procedure, one obtains T CH annel = 183 °C (shown dotted); which 
agrees favorably with the average of the measurement. Notice our calculation only provides an average temperature 
on a given surface. 

In general the channel temperature may be expressed as 

Tchannel = T A mb + Rth(T) TsupplyO —PAE) (14) 

where T AMB is the ambient (heat sink) temperature, R-m(T) is the temperature dependent total thermal resistance of 
the stack, and PAE is the power added efficiency, which may be expressed as 

„ , „ rf power developed 
PAE = — — 

^SUPPLY 

where /"supply is the DC bias power, and G p is the large signal power gain. Figure 8 gives the channel temperature 
for several cases of developed power density and values for PAE. Table 1 gives the thermal resistance breakdown 
for the layers that constitute the 8 gate device in figure 1 . Observe that the highest thermal resistance is developed in 
the GaN layer. 


1 -- 


(15) 
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Table 1 . — The change in channel temperature as the developed RF power and PAE are varied. 
This is for an 8 gate device, Lg = 0.15 pm, Wg = 100 pm, gate centered in the source-drain gap 
of 2 pm. RF powerxeff. = P D i SS 
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Figure 9. — Comparison of computed T C h a nnel with 
measured peak and simulated temperature profile 
of Park [25], versus transverse position for a 2 gate 
device for two levels of input dissipated power. 


D. Two-Finger GaN Device on Sapphire 

Measurements using nematic liquid crystal thermography on a two-gate device on sapphire were performed in 
[25]; see figure 9. Our computations yielded Tchannkl. = 59 and 76.5 °C for the two input power dissipation levels 
7 J i)iss, respectively. Again the predicted values are near a rough visual average of the channel temperature. 

Our modeling procedure can estimate the channel temperature to within 5 to 10 percent, and in some cases, 
much better than that. 


IV. SSPA THERMAL MANAGEMENT ESTIMATE 

Figure 10 is a hypothetical 50 W SSPA which consists of four MMIC chips, each of which develops 12.5 W, 
with 50 percent PAE, and they are combined by a lossless circuit. Each chip has 8 multi-gate amplifiers and each 
amplifier has 8 gate fingers. This implies the developed power density is 1.95 W/mm; which is a very reasonable 
and attainable value. The details of an amplifier on the MMIC are shown in figure 1 1(a). On the top surface the 
devices are about 360 pm long and separated 410 pm. At the heat sink their thermal footprints are about 
648x388 pm. The approximate heat spreading angles are shown dotted in figure 11(a), and the figure is not to scale. 
The thermal resistance stack is shown in figure 1 1(b), and the temperatures at the interfaces are given. Figure 11(c) 
is a cross-sectional view to demonstrate the approximate heat paths. The calculation proceeds as follows. In the GaN 
layer, the heat sources are sufficiently separated (~40 pm) so they may be treated as isolated regions. Thus the 
spreading angle is 45° andR^is 63.4 K/W using equation (7), or 56.4 K/W using equation (11). We assume the 
larger value is most accurate. The footprint at the GaN/AIN interface is 1.15x101 pm and we estimate the spreading 
angle from [26] 


tan 9 ( - k ; . 

tan 9 ; . k ; - 


(16) 
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where the angles are measured from the normal, and the subscripts i and r refer to the upper and lower layers 
respectively. This yields 9 r = 3.6°, so we can use the simple formula for R TH , that isR m = 1 /(Ak) where 1 is the layer 
thickness and ,4 is the cross-sectional area. We obtain Rjh = 12.9 K/W. The footprint at the AlN/SiC boundary is 
1.152x101 pm and we assume the 40 pm separation for simplicity. This gives Rth = 6.96 K/W and spreading angle 
of 64.8°. This implies the fluxes intersect at about 10 pm into the SiC layer. Assuming the final footprint at the 
AuSn-heatsink interface to be about 388x648 pm as in figure 1 1(c), the composite spreading angle is about 59°. 

The thermal cross-talk between adjacent regions is approximated using either [27] or [28], which yield A 7j 2 = 

6.6 or 7.4 K, respectively. Thus we choose the sink to be at 67 °C. The temperatures at the interfaces are found next. 
Starting at 7]: 7j = (1.5625 W)(0.937 K/W) + 340 K = 341. 5K. The change in effective temperature across the SiC 
layer is (1.5625)(6.96 K/W) = 10.9 K, thus the effective temperature on the top surface of the SiC layer is T = 

10.9 K + 341.5 K = 352 K. Then T A is 352.2 K using equation (15). Then T 3 = (1.5625/8)(12.9 K/W) = 354.7 K. 

The temperature change across the GaN layer is AT = (1.5625/8)(63.4 K/W) = 12.4 K. Then the effective 
temperature on the top surface is T = 354.7 K + 12.4 K = 367. 1 K. Finally, the channel temperature is obtained 
as 367.4 K from equation (14). This value, 94.4 °C is quite reasonable, and close to similar ones reported for 
similar devices. 


/ \ 



Figure 1 0. — Concept of a GaN MMIC chip with 8 amplifiers, and each amplifier con- 
sisting of 8 gates, each 100 pm wide. The four chips are assumed to be soldered 
onto a Cu-tungsten heatsink, and combined to produce nearly 50 W. 
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Figure 1 1 . — (a) Cross-section and the approximate spreading angles into each layer, (b) Corresponding thermal 
resistance stack with the interface temperatures, (c) Effective footprints at the AuSn/Cu-tungsten heatsink 
interface, which shows the 59° effective spreading angle for the entire (composite) chip. 


V. CONCLUSIONS 

We have demonstrated a simple analytical method to estimate the channel temperature in an AlGaN/GaN 
HEMT for microwave power amplifiers. It is applicable for multi-gate, multi-layer structures on SiC or Sapphire 
substrates, and for wide ranges of dissipated power. With accurate expressions for the temperature dependent 
thermal conductivity and using previously published equations for the thermal resistance of parallel stripes, one 
calculates a temperature that is within 5 to 10 percent of the true average channel temperature. 
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